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C\ ; Abstract 

The NDA C coupling constant is evaluated in a full QCD sum rule calculation. 
We study the Borel sum rule for the three point function of one pseudoscalar 
one nucleon and one A c current up to order seven in the operator product 
expansion. The Borel transform is performed with respect to the nucleon and 
A c momenta, which are taken to be equal, whereas the momentum q 2 of the 

■ pseudoscalar vertex is taken to be zero. This coupling constant is relevant in 

■ the meson cloud description of the nucleon which has been recently used to 
explain exotic events observed by the HI and ZEUS Collaborations at HERA. 
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Many explanations have been proposed for the recent observation by HI and ZEUS 
Collaborations at HERA of an excess of events at large x and Q 2 in deep inelastic scattering 
(DIS). Whereas most explanations have concentrated on physics beyond the standard model, 
some of them were based on more conventional mechanisms. One of these was suggested by 
Melnitchouk and Thomas |] • They pointed out that an enhanced charm quark distribution 
in the nucleon sea at large x could account for the observed effect. This hard component 
of the c distribution would have to be generated non-perturbatively, since the usual charm 
quark distribution arising from gluon bremsstrahlung alone looks like a typical soft sea 
distribution, namely peaks at x — > and is negligible beyond x ~ 0.4. In the works pm it 
was suggested that the charmed sea may arise from the quantum fluctuation of the nucleon 
to a virtual D~ A c configuration. A natural prediction of this picture are non-symmetric c 
and c distributions, the latter being much harder. Following this approach Melnitchouk and 
Thomas studied the effects of this charmed meson cloud on the DIS cross sections measured 
at HERA. 
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In the framework of the standard meson cloud model the anti-charm quark distribution 
in the nucleon, cjv, is given by the convolution of the valence c momentum distribution in 
the D~ meson, Cd(x), with the momentum distribution of this meson in the nucleon fo{y) '■ 

c N (x) = f dyf D (y)-c D (-). (1) 

J* y y 

The virtual D~ meson distribution in the nucleon cloud which characterizes its probability 
of carrying a fraction y of the nucleon momentum in the infinite momentum frame is given 

by I 

16tH Jo y(l-y)(s DAc - M^y \ l-y J 

where sda c = {k\ + m 2 D )/y + {k\ + M|J/(1 — y). In the above equation F(k 2 L , y) is a form 
factor at the DNA C vertex [p]]. 

In the heavy quark effective theory 0] it is assumed that the heavy quark interacting with 
light constituents, inside a hadronic bound state, exchanges momenta much smaller than 
its mass. Therefore, to a good approximation, the heavy quark moves with the velocity 
of the charmed hadron. Being almost on-shell, the heavy quark carries almost the whole 
momentum of the hadron. These considerations suggest that the c distribution in the D~ 
meson is expected to be quite hard and may be approximated by a delta function: 

c D (-)~S(--l) = xS(x-y). (3) 

y y 

From eqs. (JTJ) , (0) and (|3]) we obtain 

cvO) = Jd{x) (4) 

Similar expressions can be found for the charm quark distribution in the nucleon, cn{x). 
The above equation shows that the charm quark distributions in the nucleon are directly 
proportional to g 2 NDAc . The value used in ref. was = = -3.795. 

Plugging the above calculated cn (x), toghether with parton distributions of all flavors, 
in the standard DIS cross sections expressions, da j dxdQ 2 , Melnitchouk and Thomas were 
able to reproduce the observed rise of the cross sections at large values of M (M ~ 200 GeV) 
and Q 2 . 

In this work we shall concentrate our attention on the coupling constant gNDA c , since it 
enters directly in the calculation of the DIS cross section. The value used to reproduce the 
HERA events in |l| , while a reasonable guess, has neither a solid theoretical justification 
nor a phenomenological basis. It may, however, be calculated in the framework of QCD 
Sum Rules (QCDSR). In the literature there are many calculations of coupling constants 
with QCDSR. We shall make reference to two of them: to |J, which describes the general 
method and to 0, which is recent and estimates gNKK- 

In order to calculate the N D~ A c coupling constant using the QCD sum rule approach 
we consider the three-point function 

A(p,p',q) = ! d^xd 4 y(0\T{ VA Xx) j5 (y)rJ N (0)}\0)e^e-^ , (5) 
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constructed with two baryon currents, t]a c and 7] n, for A c and the nucleon respectively, and 
the pseudoscalar meson D current, j'5, given by |||| 

Vk c = £dbc(UaC f y 5 d b )Q c , (6) 

VN = £abc(Ua C 1^)151^0 , (7) 

35 = Qijsu , (8) 

where Q, u and d are the charm, up and down quark fields respectively, C is the charge 
conjugation matrix. 

Due to restrictions from Lorentz, parity and charge conjugation invariance the general 
expression for A(p,p',q) in Eq.(|5|) has the form 

A(p,p', q) = Fi(p 2 ,p' 2 , g 2 )27 5 + F 2 (p 2 ,p 2 } q 2 )^ 5 



+ F 3 (p 2 ,p' 2 , q 2 )Pz l5 + F A (p 2 ,p' 2 , q 2 )a^ 5 p,p' v , (9) 

where q = p' — p and P = (p + p')/2. 

In the phenomenological side the different Lorentz structures appearing in Eq.@ are 
obtained by the consideration of the A c and N states contribution to the matrix element in 
Eq. (|): 

(0|r /Ac |A c (pO)(A c (p / )|j5|iV(p))(iV(p)|^|0) , (10) 
where the matrix element of the pseudoscalar current defines the pseudoscalar form-factor 

(A c (p')\j 5 \N(p)) = g P {q 2 )u{p')i^u{p) , (11) 
where u{p) is a Dirac spinor and gp(q 2 ) is related to gNDk c through the relation |j| 

D 



(m u + m c ) q 2 — m 



where and fo are the meson D mass and decay constant and m u ( c \ is the u(c) quark 
mass. 

The other matrix elements contained in Eq.([10|) are of the form 

(o\ Va Mp')) = *aXp') (13) 
(N(p)\rj N \0) = \ N u(p) , (14) 

where Aa c and A at are the couplings of the currents with the respective hadronic states. 

Saturating the correlation function Eq.(f|) with A c and N intermediate states, and using 
Eqs.flTTj), (0), (0), (|TJ) and fllj we get 

^'■«> = ^ ^lt^i f /-K ^^ + hl8te " ■ (15) 
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which can be rewritten as 



A^ hen \p,p',q)=\ Ac \ 



2m 2 D f D 9nd\ c 



+ 



(m u + m c ) q 2 — m 2 D p' 2 — M A p 2 
M Ac + M N 



M 2 N 



[(M Ac M N - p.p')t l5 



-#7b + (M Ac - M N )Pt l5 - o» v lhVilV ' v 



+ higher resonances , (16) 



where we clearly see all the Lorentz structures present in Eq.([|). We will follow refs. |||| 
and write a sum rule for the structure $75. As we are interested in the value of the coupling 
constant at q 2 = 0, we will make a Borel transform to both p 2 = p' 2 — > M 2 . The contribution 
from excited baryons will be taken into account as usual through the standard form of ref. 
0. 

As in any QCD sum rule calculation, our goal is to make a match between the two 
representations of the correlation function (J^) at a certain region of M 2 : the OPE side and 
the phenomenological side. The contributions to the OPE side are represented in figures 1, 
2 and 3. 

In the OPE side only odd dimension operators contribute to the ^75 structure, since the 
dimension of Eq.(|5|) is four and $ takes away one dimension. Therefore, the perturbative 
diagram in Fig. 1 contributes through the m c operator. To evaluate the perturbative 
contribution we write a double dispersion relation to the amplitude F2 in Eq.(|9|) and use the 
Cutkosky's rules to evaluate the double discontinuity (see ref. 0]). After Borel transforming 
with respect to P 2 



-p = —p' , and subtracting the continuum contribution, we get 



F 2 (M 2 ,q 2 ' 



1 ruo 
An z Jml JO 



dsp(u, s, q 2 



( e -u/M 2 _ e -s/M 2 ^ 



S — U 



(17) 



with 



p(u,s,q 2 



3 m r 



2(2vr) 2 ^X(s lUl q 2 ) J ° 



dm <m — h 



4 2A(s,w,g 2 ) 



m 
T 



1 

— + 



y/s, 



&((cos6 K y - 1) 



2s \(s,u,q 2 ) \ 

where p' Q = (s + u — q 2 ) / (2 y/s) and 

u + m 2 — m 2 — p' Q (s + m 2 )/ y/s 



(18) 



cos 9 



K 



2s- 



(s — m 2 )J\(s } u, q 2 ) 



(19) 



In Eq. (pT|) i*2 stands for the Borel transformation of the amplitude F2, the subscript 1 refers 
to the diagram 1 and uq gives the continuum threshold for A c . 

The next lowest dimension operator is the quark condensate with dimension three (Fig. 
2). In ref. || only the diagram shown in Fig. 2a was considered because the pion mass 
was neglected and, therefore, only terms proportional to 1/q 2 contributed. Since in our case 
mjj will not be neglected we have also to consider the diagrams in Figs. 2b and 2c. Their 
contributions are given by 
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F 2 (M\q*) 



(qq) 



2a 



?7r 2 q 2 — m 2 



M A E 1 



(20) 



with Ei = l- e - u °l M \l + u /M 2 ). 

(qq) r«o 



26+2c 87T 2 

where 5 = 1/2 — a(u + s)/(2s) and 



u " mc ^ (2-<5 + ft)s + (l + o> ^_ M/M2 _ /M 2 
( M -s) 2 16 



, (21) 



2s 



[u - s 



m r + 



s — u 



(22) 



The last class of diagrams which we will consider is the dimension 7 operators of the type 
m c (qqqq) (q being &u or d quark) shown in Fig. 3. Other dimension 7 operators come from 
graphs which contain at least one loop and are suppressed by factors 1/47T 2 . The expressions 
for these contributions are 



F 2 (M 2 ,q 2 " 



3a+3b 



(qqqq) m c 
6 q 2 — m 2 c 



(23) 



F 2 (M 2 ,q 2 " 



(qqqq) 1 - e" m '/ M2 
3c 6 m c 

The Borel transformation of the phenomenological side gives 

2 )JD yNDAc 

m u + m c ) q 2 — m 2 D M\ c — M 



F 2 (M 2 ,q 2 ' 



\ x mD ^ D 9nda c M Ac + M N M 2 r/M 2 
(m„ + ran) a 2 — m 2 r> Ml — M% 



M 2 JM 2 ^ 



e A <= 



(24) 



(25) 



For Aa c and A a? we use the values obtained from the respective mass sum rules for the 
nucleon [|]|| and for A c Q: 



I x I2.-M2 /M 2 



ml 



5127T 4 



-u/M 2 



ml 



u u 



m„ 



1 £ 1 + — - 12 In — £ 



m 2 ' mf. 



+ 



-m 2 JM 2 



(26) 



IXn1 2 6 -M 2 n /M 2 



M 6 n 2 /- - V 

^ 2 + -(qqqQ) 



(27) 



where E 2 — 1 — e _s °/ M2 (l + Sq/M 2 + Sq/(2M 4 )) accounts for the continuum contribution 
with So being the continuum threshold for the nucleon. We have neglected the contribution 
of the gluon condensate in the mass sum rules and in the three point function since it is of 
little influence. 



In order to obtain gNDA c we identify eq. ( p5D with the sum of eqs. (jl7|, ^T], |23], [24]) and 
using eqs. fl2"E| , f2"7| ) we solve the resulting equation for gNDA c as a function of the Borel mass 
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squared. The resulting curves are shown in Figs. 4 and 5. Since eqs. ([261 , W^i ) determine 
only the absolute value of and Aa c we can not determine the sign of Qnda c - 

In this calculation there is no free parameter. There are, instead, some numbers which are 
heavily constrained by other theoretical or phenomenological analyses. They are the sources 
of uncertainties in our result. The quark condensate was taken to be (qq) = — (0.23) 3 GeV 3 . 
The continuum thresholds appearing in E\ and Ei were chosen to be respectively uq = 
(M Ac +0.7) 2 GeV 2 and s = (M N + 0.7) 2 GeV 2 . The hadron masses are M N = 0.938 GeV, 
Ma c = 2.285 GeV and vtld = 1.8 GeV. The charm quark mass was taken to be m c = 
1.3 and 1.5 GeV. In order to take into account possible deviations from the factorization 



hypothesis []T2| we introduce a correction factor K in the formula 

(qqqq) = K(qq) 2 (28) 

and consider the cases K = 1 and K = 2. Finally we need the D meson decay constant 
fo- This constant will be probably evaluated in the next generation of experiments on fully 
leptonic D decays. So far it has been usually estimated |13|] to be fjj = 200 ± 30 MeV. 



On the other hand, new measurements of fo a [1^] by the CLEO Collaboration point to 
f Ds = 344 ± 37 MeV. Considering that f D and f Da are related by /dJ/d ~ 1.15 - 1.2 
this indicates that fjj ~ 300 MeV. In view of this we take fjj in the interval 200 — 300 MeV. 

The relevant Borel mass here is M ~ m n+^Ma c an( j th ere f ore we analyse the sum rule 
in the interval 1.5 < M 2 < 3.5 GeV 2 . In Fig. 4 we plot gNDA c as a function of M 2 for 
different values of fjj and K. As it can be seen, the curve is rather flat showing that, for 
QCD sum rules standards, this sum rule result is stable. Moreover, the uncertainties of 50% 
both in f D and K lead to a spread in the value of gNDA c - In Fig- 5 we fix f D = 250 MeV 
and K = 1 and vary m c from 1.3 (solid line) to 1.5 GeV (dashed line). The corresponding 
variation in gNDA c is of the order of 15%. We have also varied the continuum thresholds 
in the interval M N + 0.5 < s < M w + 1.0 GeV and M Ac + 0.5 < u < M Ac + 1.0 GeV. 
Under these variations the final result does not change appreciably. This indicates that the 
continuum contribution is under control. 

Considering all the uncertainties discussed above our final result is 



\9nda c 



47T 



1.9 ± 0.6 (29) 



This value is smaller than the one used in the calculations of ref. [0 by a factor two. 
At first sight, since the normalization of the c distribution in the nucleon is proportional 
to the square of gNDA c , using our value would significantly reduce the contribution of the 
non-perturbative nucleon charm sea employed to study the HERA events. However there 
are many other possible intermediate states apart from D — A c which contribute to generate 
this charm sea like, for example, D* — A c , D — S and D* — E. These higher mass states are 
neglected in some calculations but, as shown in ref. [[|, already in the strange sector they 
may give important contributions for the evaluation of strange form factors and strange 
quark momentum distributions. In the charm sector we expect the contribution of these 
states to be even more important. Our conclusion is therefore, that the non-perturbatively 
generated charm sea remains qualitatively a good candidate to explain the observed exotic 
events. Quantitative calculations should, however, include the contribution of other charm 
mesonic clouds, specially in view of the small value found here for gNDA c - 
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Figure Captions 



Fig. 1 Perturbative contribution to the OPE. 

Fig. 2 Diagrams that contribute to the Wilson coefficient of the operator qq. 

Fig. 3 Diagrams that contribute to the Wilson coefficient of the operator qqqq. 

Fig. 4 Qnda c as a function of the squared Borel mass M 2 . The different lines correspond 
to different choices for the D meson decay constant, /d, and for the factorization 
correction factor K as indicated in the legends. 

Fig. 5 The same as Fig. 4 for diferent values of the c quark mass, fo = 250 MeV and K — 1. 
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fig. 2 
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